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Manipulation of skyrmions in nanodisks with a current pulse and skyrmion rectifier 
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A skyrmion in a nanosized disk of a chiral magnet can be used as a bit of information. To this end, it is 
desirable to control the creation and removal of a skyrmion only by currents without using external magnetic 
fields. Here we propose to create a skyrmion by applying a current pulse to a nanodisk. The skyrmion can be 
removed from the disk by applying a dc current. We show that the dynamics of the created skyrmion can lead 
to a rectification effect, in which a dc voltage is generated by the motion of skyrmion in the presence of an ac 
current. 

PACS numbers: 75.10.Hk, 75.25.-j, 75.30.Kz, 72.25.-b 



A new spin texture called skyrmion has been found in cer- 
tain chiral magnets without inversion symmetry, such as MnSi 
and Feo.5Coo.5Si. fl]-0 A skyrmion can be described in terms 
of spins which wrap a sphere. When a conduction electron 
passes through the skyrmion, its spin is fully polarized by the 
spin texture of the skyrmion, yielding a quantized Berry phase 
An that corresponds to a quantized magnetic flux O = hc/e. 
This emergent electromagnetic field significantly affects the 
dynamics of the skyrmions and produces a Magnus force per- 
pendicular to the skyrmion velocity. (8J [9) Application of 
a spin polarized current exerts a spin transfer torque on the 
skyrmion, causing it to move at an angle to the current direc- 
tion determined by the dissipation. ifTUl The resulting voltage 
is perpendicular to the velocity; thus, a Hall (transverse) volt- 
age is induced, which has been measured experimentally. IfTTI 

Since skyrmions are stable topological excitations in chiral 
magnets, they have promising applications in spintronics. Ifl2l 
Due to the presence of the Magnus force (DO, the thresh- 
old current to drive the skyrmion from the pinning sites is ex- 
tremely small, 4-5 orders of magnitude smaller than that for 
magnetic domain walls. IfTTI [T3l [T4ll Thus, skyrmions have 
tremendous advantages for information storage applications, 
because the lower current implies less Joule heating. From a 
technological point of view, it is desirable to be able to ma- 
nipulate, create, or destroy skyrmions by applying only a cur- 
rent without using an external magnetic field. This requires 
use of magnetic disk that is is typically the same size as the 
skyrmion, of the order of 100 nm, such that only one skyrmion 
can be accommodated in the disk. The goal of producing a 
single magnetic vortex in a nanodisk was recently achieved 
|fT5l[T61 . In most known chiral magnets, the skyrmion is not 
stable without an external magnetic field, and in the absence 
of a field a spiral structure is favored. Upon increasing the 
magnetic field, a skyrmion state becomes the ground state. 
[ 17 1 Two immediate questions regarding the dynamic creation 
of skyrmions arise: first, how to stabilize skyrmions without 
magnetic fields, and secondly, how to create a skyrmion using 
currents. 

The answer to the first question is to exploit the magnetic 
anisotropy in thin films. 1 18 1 Due to the reduced dimensional- 
ity, an out-of-plane anisotropy develops that is inversely pro- 
portional to the film thickness. Such anisotropy to some ex- 



tent acts as a perpendicular external field, so that skyrmions 
are easily stabilized in thin films. |5][T5-22| For the second 
question, it was demonstrated in Ref. 23 that skyrmions can 
be created by applying a circular current, which acts as an 
effective local magnetic field. However, implementation of 
a circular current may be difficult in many practical applica- 
tions, especially for disk sizes that are below one micrometer. 
Here we introduce a mechanism of destabilizing the ferromag- 
netic (FM) state with a transport current, and then creating a 
metastable skyrmion. 

We also explore the dynamics of the created skyrmion in 
a nanodisk. The Magnus force F,„ ~ zx vis always per- 
pendicular to the velocity, so that the effective dynamics of 
the skyrmion is equivalent to a charged particle moving in a 
transverse magnetic field. Here I is a unit vector perpendicu- 
lar to the disk. The trajectory of the skyrmion is thus chiral, 
and is either clockwise or counterclockwise depending on the 
sign of the Magnus force. Under an ac current, we show that 
the trajectory of the skyrmion is an ellipse in the presence of 
a circularly symmetric confining potential. Since the voltage 
is perpendicular to the velocity, one thus can use skyrmions to 
extract a dc voltage by applying an ac current, indicating that 
it is possible to create a skyrmion rectifier. 

We consider a thin circular disk of chiral magnet with typ- 
ical size of the order of 100 nm, comparable to the size of a 
single skyrmion. The system is described by the Hamiltonian 
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where the first term is the exchange interaction, the second 
term is the Dzyaloshinskii-Moriya (DM) interaction due to the 
spin-orbit coupling [28-30], and the third term is the magnetic 
anisotropy with the easy axis perpendicular to the disk. For a 
thin film of thickness d, we have J a °c l/d. [181 At low tem- 
peratures, the amplitude of the magnetization M is approxi- 
mately fixed and is normalized in units of n = M/M s with M s 
the saturation magnetization. The dynamics of the spins are 
governed by the Landau-Lifshitz-Gilbert equation of motion 
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caused by an applied spin polarized current. [36] The spec- 
trum of the spin wave excitation in the presence of a dc current 
J and the dissipation is 



FIG. 1. (color online) (a) Schematic view of a chiral magnet sub- 
jected to a current. To study the skyrmion rectification effect, we 
measure the voltage along the lines shown in the figure when an ac 
current is injected, (b) A twisted magnetic spiral at J A = 0.4J ex /D 
and (c) a metastable skyrmion at J A = 3.4J ex /D. The vectors in 
the plots denote the n x and n y components, and color denotes the n z 
component. Here the radius is R = S.OJ^/D. 



where the effective local field is H ei f = -blrilbw - 7 ex V 2 n - 
2DV x n + Jah z z. The first term is the spin transfer torque and 
the last term is the Gilbert damping with coefficient a. Here 
y - a 3 litis) with a the lattice constant and s the total spin. The 
ground state can be found by numerical annealing by adding a 
Gaussian noise field to H e g. For weak anisotropy, the ground 
state is a magnetic spiral, as shown in Fig. [TTb). In the con- 
fined geometry, the direction of the spiral is twisted. Upon in- 
creasing the anisotropy, a FM state is stabilized. In this system 
without an external field, the out-of-plane anisotropy does not 
stabilize a skyrmion in the ground state. The skyrmion can, 
however, be a metastable state. We initially put a skyrmion at 
the center of the disk, and find in our simulation [34| that the 
skyrmion is stabilized for 2.5D 2 /J ex < J a ^ 4.5D 2 /J ex . For 
a large J a, the size of the skyrmion shrinks because the whirl 
of spins would cost more energy in this case. The topological 
charge for the skyrmion is Q — J dr 2 q(r) with the skyrmion 
density q - n-(3 x nx5 v n)/(47r). [35 1 The electric field induced 
by the motion of the skyrmion is E = tax ■ (Vn x d,n) j2e. ifTOl 

We next study the creation of a skyrmion by a current pulse. 
The following conditions must be met. First, the anisotropy 
J A should be in the range 2.5D 2 /J ex < J A < 4.5D 2 /J ex to 
support a metastable skyrmion excitation. Secondly, the am- 
plitude of the current must be large enough to trigger the in- 
stability of FM state, as detailed below. Thirdly, the duration 
of the pulse should be large enough to allow the system to 
evolve into the skyrmion state. Typically this time is of the or- 
der of !QJ ex /(yD 2 ). Finally, the damping constant should be 
moderately large. In simulations a skyrmion can be created 
easily with a > 0.1. For small damping, it is more difficult 
to stabilize the created skyrmion inside the disk because the 
skyrmion escapes from the disk during the time evolution. 

It is possible to create a skyrmion using current J because 
the FM state can be destabilized by a spin transfer torque 
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The imaginary part of Im[Q] is due to the Gilbert damp- 
ing. The spectrum becomes gapless when J is larger than 
J > Ji - 2yypA~Ja./(a 2 + 1), which indicates an instabil- 
ity of the FM state. In the presence of weak perturbations, 
such as thermal/quantum fluctuations of spins, the perturba- 
tions increase exponentially. The external current couples to 
the emergent vector potential A associated with spin n, where 
A = ichb^Vb/e with b the spin coherent state [35|. Thus the 
state with a non-zero vector potential A is favored after the 
ferromagnetic state has been made unstable. In the presence 
of the DM interaction, the state with a non-zero emergent vec- 
tor potential is the skyrmion state. 

To demonstrate the dynamical creation of skyrmions by ap- 
plying a current pulse, we perform numerical simulations of 
Eqs. ([T) and (|2]i by adding an arbitrarily weak white noise 
magnetic field to H e ff. For numerical details, see Ref. 11341 . 
During the time evolution when J > /,■ as shown in Fig. [2j 
the x or y component of the spin increases, and a skyrmion is 
created in the final state. The created skyrmion is then driven 
by the current, moves towards the edge of the disk, and finally 
disappears. This processes repeats if a dc current is applied. 
In the presence of a current pulse, the created skyrmion will 
remain inside the disk if the current is turned off before the 
skyrmion has time to leave the disk. 

The nanodisk provides a geometric confinement that orig- 
inates from the boundary condition, which prevents it from 
leaving the disk when no current is injected. This geomet- 
ric confinement is effective only when the skyrmion is near 
the edge of the disk. In simulations, we can introduce a pin- 
ning center by modulating the electronic density at the center 
of the disk, since the exchange interaction is proportional to 
the electronic density. When J ex {r) is depressed in the center 
of the disk, it provides a pinning potential for skyrmion. [9] 
In the case when the center of the skyrmion is already out- 
side the disk, the skyrmion will be pushed away from the disk 
even though the current is turned off. In this situation, one 
may repeat the application of a current pulse with a slightly 
different duration to create a skyrmion. The removal of the 
skyrmion can be achieved by applying a dc current with am- 
plitude smaller than /, but larger than the depinning current 
corresponding to the pinning potential. In this way, one can 
create and remove a skyrmion in a nanodisk in a controlled 
way using only a current. 

In the absence of an external magnetic field, the states ob- 
tained by reversing all the spins along the easy axis (the z 
axis here) are degenerate. Therefore there are two types of 
skyrmions. In one type of skyrmion configuration, the spins 
in the skyrmion center are down and the spins away from the 
center are up (type I). The other degenerate skyrmion has the 
reverse spin direction, i.e., spins in the skyrmion center are 
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FIG. 2. (color online) (a) Profile of the current pulse used in simulations, (b)-(f) Time evolution of the spin texture in the disk. The vectors in 
the plots denote the n x and n y components, and color denotes the n- component, (b) and (c): After applying the current pulse, an instability in 
the FM state is triggered, (d): Three skyrmions are created after turning off the pulse, with the center of one skyrmion outside the disk, (e): 
The skyrmion at the edge in (c) leaves the disk and the other skyrmion is repelled towards the edges due to the inter- skyrmion repulsion, (f): 
Finally, one skyrmion is stabilized at the center of the disk. The radius of the disk is R = &.0J cx /D and J A = 3.0D 2 /J SX . The time in the plot is 
in units of t = J ex /(yD 2 ) and current is in units of J = 2De/Ti. For MnSi, t * 0.02 ns and J ss 10 12 A/m 2 . To pin the created skyrmion, we 
introduce a defect at the center of the disk by assuming 7 cx (r) = / ex [l - 0.8 exp(-r)]. In simulations, an arbitrarily weak Gaussian noise field 
is added to H c jf to trigger the instability. 



up and the spins away from the center are down (type II). 
However, the chirality for both skyrmion configurations is the 
same, since it is determined by the sign of the DM interaction. 
In simulations, we found that if one starts from the FM state 
with spin up, one creates a skyrmion of type I, and vice versa. 
Here we investigate the rectification effect, or the conver- 
sion of an ac current into a dc voltage, that originates from 
the unique dynamics of skyrmions in a nanodisk. To demon- 
strate the idea, we consider the particle-level description of a 
skyrmion in the presence of a pinning potential in dimension- 
less form |9| 



77v = |xv + zxJ + ¥ d (r < - r,), 
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where the first term on the right-hand side is the Magnus force, 
the second term is the Lorentz force, and the last term is the 
pinning force Fj - — VC7(r ; - - r,) with a pinning potential 
U(Tj - r,). Here the damping coefficient 77 has the contri- 
butions from the Gilbert damping and the dissipation due to 
the conduction electrons. The unique feature of the dynamic 
equation is the presence of the Magnus force, which causes 
the skyrmion motion to be clockwise in the x — y plane. The 
trajectory of a skyrmion moving in a strong pinning potential 
and subject to an ac current J x - 7 x o sin cot forms a limit cycle, 
according to the Poincare-Bendixson theorem, |37] because 
(1) there is no fixed point; and (2) the skyrmion trajectory 
is bounded by the pinning potential. For a parabolic pinning 
potential U - U p r 2 /2, the trajectory of the skyrmion can be 
found analytically and is an ellipse: 



skyrmion velocity in the azimuthal direction V r oc v^, one can 
extract a dc voltage even though only an ac current is applied, 
i.e. a skyrmion rectifier appears. 

We consider a nanodisk described by the Hamiltonian in 
Eq. ([TJ, which supports a skyrmion as a metastable state. We 
explicitly introduce a defect at the center of disk to pin the 
skyrmion by assuming J ex (r) = J ex [l _ 0.8 exp(-r)]. Initially 
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When Up ~ co and r\ <k 1, the trajectory is nearly circular, 
and when 77 » 1, it is nearly linear. If one measures the volt- 
age along the radial direction V r , which is proportional to the 



FIG. 3. (color online) Dependence on a of (a) the radial voltage 
V x measured for < X < R at y = and (b) the radial voltage V y 
measured for 0<y<^atx = 0at different driving frequencies 
a> of the ac current, (c) and (d) The corresponding trajectories of 
skyrmions at different a and a>. An ac current along the x direction 
is injected with the form J x = J x0 sm(<i)t), where the amplitude J M = 
0.2De/%. A defect is introduced at the center of the disk by assuming 
/«(r) = /„[!- 0.8 exp(-r)]. 



we put a skyrmion in the center of the disk and add an ac cur- 
rent with amplitude smaller than the depinning current, such 
that the skyrmion is not driven outside the disk. We calculate 
the time averaged voltage between {R, 0) and (0, 0), as well as 
the voltage between (0, R) and (0, 0), as shown in Fig. [T] by 
solving Eqs. ([TJ and (|2]) numerically Il34ll . We also calculate 
the time dependence of the center of mass defined as 



hi 



dr q(r)r. 
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The trajectories at several different values of a are shown in 
Fig. [3] As a increases the dissipative force becomes dominant 
over the Mangus force, causing the trajectory to become ellip- 
tical. In the a — > oo limit, the trajectory becomes a straight 
line in the radial direction and thus the time averaged dc volt- 
age vanishes. The dependence of the voltage on a is shown in 
Fig. [3] There is a clear correspondence between the voltage 
and the skyrmion trajectory. For instance, as a increases, the 
larger principal axis of the ellipse approaches they axis. Thus 
the amplitude of V t (0 < x < R) increases while V v (0 < y < R) 
decreases in most cases. Thus the voltage in Figs. [3] (a) and 
(b) can be inferred from the trajectories in Figs. [31(c) and (d). 
Particularly for a> - 0.1, the dependence of V x on a is non- 
monotonic. All these behaviors are captured by the particle 
model in Eq. Q, 

We use typical parameters for MnSi to estimate the opti- 
mal parameters for the creation of a skyrmion. For MnSi, 
/ex ~ 3 meV/a and D » 0.3 meV/a 2 with the lattice constant 
a « 2.9 A. 1 10] For a disk size of about 100 nm, the strength 
of the current pulse should be larger than J > /,■ ~ 10 12 A/m 2 , 
while the duration should be longer than 0.2 ns. To reach the 
range of the uniaxial anisotropy 2.5D 2 /J ex < Ja < 4-.5D 2 /J ex 
required for stabilizing a skyrmion, the thickness of the disk 
should be around d « 5 nm. Ifl8l One challenge here is to en- 
hance the damping coefficient a > 0. 1 for a skyrmion, since 
the bare Gilbert damping is usually very weak, ac < 0. 1 . Nev- 
ertheless, there is an additional damping mechanism due to the 
electric field induced by the motion of a skyrmion. 1 10 1 The 
dressed damping thus can be larger than a > 0.1, which is 
optimal for the creation of skyrmions. 

Finally, we estimate the conversion efficiency of the 
skyrmion rectifier, which is defined as the ratio of the induced 
dc voltage Vd c to the amplitude of the ac current P = Vdc/Ao- 
According to the results in Fig. [3] for a single skyrmion, we 
obtain P w 10~ 17 V ■ m 2 /A. The rectification effect is weak; 
however, it is experimentally observable. For an applied cur- 
rent 7 x0 ~ 10 12 A/m 2 , we have V dc w 10 //V. Note that V dc 
can be enhanced at a given 7 x o by using many skyrmions in a 
large disk, since the voltage is proportional to the number of 
skyrmions. In a large thin film without pinning centers, Vh c is 
just the Hall voltage associated with the motion of skyrmions, 
which has been measured experimentally. ifTTll 

The applied current generates a magnetic field and we esti- 
mate its strength. For a current density 10 12 A/m 2 , disk radius 
100 nm and disk thickness 5 nm, the induced field is about 
25 G at the surface of the disk. The generated field is much 



weaker than that used in experiments to stabilize the skyrmion 
phase. Thus we expect the self-induced magnetic field has 
a negligible effect on the dynamical creation and removal of 
skyrmions, and also the skyrmion dynamics discussed here. 

To summarize, we have demonstrated a method to create 
and annihilate skyrmions in chiral magnetic nanodisks with- 
out an applied magnetic field by utilizing an applied current 
to create an instability in the ferromagnetic state. The abil- 
ity to stabilize skyrmions in nanodisks would be a significant 
step forward for possible applications of skyrmions. Addition- 
ally, we show that once a skyrmion is stabilized in a disk, the 
unique skyrmion dynamics produce complex skyrmion orbits 
under the application of an ac drive, leading to the production 
of a dc voltage similar to a rectification effect. 
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